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Abstract 

We present estimates of the charming penguin contribution to B — > Kir, tttt, Kn, Krj' 
decays due to intermediate charmed meson states. We find that this contribution 
is indeed significant for B — > Kir decays, and its inclusion, together with the tree 
and penguin terms, produces large branching ratios in agreement with data, though 
the analysis is affected by large theoretical uncertainties. On the other hand, for 
B — ► 7T7T, Kn, Kn' decays, the effect of the charming penguin contribution is more 
modest. We also compute CP asymmetries for B — ► Kir, tttt decays and we obtain 
rather large results. 



PACS: 13.25.Hw 



1 Introduction 



In a recent paper |TJ , hereafter referred to as X, we gave an estimate of the so-called charming 
penguin [[|, |§ contributions to the decays B — > Kn. This is a long-distance part of the 
decay amplitude whose imaginary part results from the decay chains 

B -> D(D S ) -> Ktt , 

£ -> D*(D* s )^Kir, (1) 

while the real part can be computed by a tree diagram of the effective chiral lagrangian 
for heavy mesons @]-|§. In the present paper, we shall call this amplitude AchP- The 
relevance of these contributions for B — > Ktt decays was first pointed out in JXO | ; though 



suppressed in the factorization approximation, these terms are enhanced by the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix factor V c bV* s in comparison to the short distance terms, 
i.e., contributions arising from the Tree and Penguin terms in the factorization approxima- 
tion, whose amplitude we call here At+p- In X we have shown that, even taking into account 
the uncertainties inherent to this calculation, the contribution of the charming penguins con- 
tributions to the decay channels B + — > K°n + and B° — > K + ir~ is indeed significant and can 
explain the difference between the data and the result obtained by At+p- In the present 
paper, we wish to extend the analysis to cover other B decay channels with a Ktt pair in 
the final state, as well as other charmless B decays into two pseudoscalar mesons, i.e., 

B -> 7T7T , (2) 

B -> Krj, (3) 
B -> Krf . (4) 

For the processes (|J) and (^) we add to At+p the the charming penguin contributions 
with D^*\ Dy intermediate charmed meson states; for the irir final state, the charming 
penguin contribution is obtained by D^*\ intermediate states. All the relevant formulae 
are presented in X and can be applied here with some obvious changes. For example the 
substitution K — > 7r for the channel or the substitution of the pion physical constants 
with the analogous observables of r] and rj' for the channels (|3|) and (ffl). A few points, 
however, deserve a more detailed discussion; let us examine them in the next section. 



2 Discussion on the method and its uncertainties 

The procedure for obtaining the real part is based on the use of an effective field theory 
satisfying chiral symmetry as well as heavy flavor symmetries. The main point in this 
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procedure is the following approximation (we take B + K°n + as the representative channel 
for the B — > PP decays): 

A Ch p = ^a 2 V; b V cs <K n + \:J,(0)J^0):\B + >^ 

« ^a 2 V; b V cs J^<K\ + \T{j,(x )nO)}\B+> (5) 

with Jn = 67^(1 — 75)0 and = cry v {\ — 7s)s; moreover a 2 = (c 2 + ci/3) (= 1.03) where Ci 
and C 2 are Wilson coefficients, and 

x x = (0, n//i) , (6) 

where \n\ = 1 and \i is a scale representing the onset of the scaling behaviour. In (|5|) we 
have not considered color octet operators which would have given no contribution as we 
consider only the color-singlet physical intermediate states. This approximation is based on 



the light-cone expansion [11], [O], which in the present case reads: 



< K°7T + \T { J M (x )>(0)} \B + > w < K°7T + \ : J M (0)>(0) : \B+ > + O(zg) , (7) 

where the 0(xq) terms are negligible for \i = 1/\xq\ sufficiently large (fi ~ m&); clearly 
the integral over n corresponds to an average over the directions of x$. Moreover, the non 
trivial scale dependence of the r.h.s. of eq. (^) is matched by the Wilson coefficient, to 
give scale-independent physical observables. Also the l.h.s. contains a scale (the cut-off). 
Ideally this scale dependence should be cancelled by the short distance coefficient as well; 
in practice, however, the cancellation is not complete as we make a truncation in the long 
distance physics and we include only the low lying charmed intermediate states. In order to 
compute eq. ([5]), we write 



<^°vr + |T{j,(x )>(0)}[fi + >= Jj^e-^T(Q), (8) 



where 



T(q) = J d 4 xe +iq - x < K°7T + \T{j^(x)J^0)} \B+ > . (9) 

Let us now show that, after averaging over n, one obtains a cutoff over the high frequencies 
in eq. (§). As a matter of fact, one has 

AchP = ^vsy. j$L-m jt = %^ v - !ik T(q)e( »< ' 

(10) 

where the cutoff function is 

sin|g1//i 

OQi, WD = ■ (11) 
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For fi — > oo, 9(fi, \q\) — > 1; for finite values of /i, this function cuts off from the q— integral 
in eq. ( |i~0| ) the region — q 2 > fi 2 . Instead of the smooth oscillating function (pTTj) , we used in 
X the step function 

% 2 + ^ 2 ), (12) 

which allows a considerable reduction of computing time. We also stress that one can 
extract from the integration in the momentum q the heavy mass contribution according to 
the formula 

Q = Pb ~ P D (*) = ( m B ~ m D (*))v - t (13) 

(see Eq.(30) of X). Here is the heavy meson velocity and I is a residual momentum. By 
this the cutoff function on the i— integration becomes 

9(f + A) , (14) 

and the value of the cut-off fi£ found in X is 0.5 — 0.7 GeV. The whole procedure we have 
described so far has been used several times in the past in the application of the light cone 
expansion ideas to the nonleptonic weak decays, starting from the pioneering work of K. 



Wilson (see Section 7 of |LT| and the subsequent work of several authors |[T3[). We repeat 
it here as it may not be familiar to some readers and also to stress that the correct value 
of the cut-off /i (corresponding to fii) is not the W mass, but a scale of the order of m h or, 
better, mj — m c . It is a consequence of the precocity of the scaling behaviour, a well-known 
example of which is provided by deep inelastic scatteringfj. The introduction of non-locality 
through the cutoff /i of the order of nib (or fi£ of the order of 0.5-0.7 GeV) has a clear physical 
meaning. It corresponds to a separation between short distance physics (whose physical fea- 
tures are embodied in the Wilson coefficients) and long-distance physics, which is dominated 
by hadronic states and resonances. In considering the long distance part, we have included 
the low-lying states that could contribute, i.e., D^*'Dy . This procedure can be avoided for 
operators having factorizable contributions; in the case of the charming penguin contribu- 
tion, however, such contributions do not exist, while the non factorizable contribution are 
Cabibbo enhanced. This is the reason for taking them into account explicitly. 

The second point to be stressed is that while in the present paper, as well as in X, we are 
using the chiral lagrangian effective theory for heavy mesons 0-0, the light pseudoscalar 
mesons in the final state have large momenta. Therefore, the effective lagrangian must be 
corrected to take into account the hard meson momenta. The procedure we adopted in X 
was based on the introduction of form factors, similar to the approach followed in 0, [0 
and |]]. We were able to estimate, by using the constituent quark model, the form factor 
correcting the B*Bir and D*Dn vertices. In addition, one should also consider the form 



*On the basis of the previous remarks, the criticism of I contained in |l4| appears to be unjustified. 
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factor correcting the diagram in fig. 1, which represents the main contribution to the real 
part of AchP ( m this diagram is depicted in fig. 2a). The weak vertices — > Kir 



(I v 




Figure 1: The charming penguin contribution diagram dominating the real part of the amplitude. 
The boxes represent weak couplings. The DKir and the D*Ktt couplings represent the direct 
coupling of the effective chiral theory for heavy mesons. 

correspond to direct, non- resonant couplings and arise from the weak effective current: 

L, a = y Tr^( 1 - 75 )H b £l a , (15) 

which is the effective realization of the quark current g a 7 M (l — 7s)Q. ct is related to the heavy 
meson leptonic decay constant by the formula a = fo ^mp, valid in the infinite quark mass 
limit. Moreover, 

Ha = l Ar- (K^ - Pall) (16) 



and 



£ = e l ~ . (17) 

In these formulae, v is the heavy meson velocity; P a , P*^ are the annihilation operators of 
heavy pseudoscalar and vector mesons made up by a heavy quark and a light antiquark of 
flavour a (a = 1,2,3 for u, d, s); Ai is the usual 3x3 matrix comprising the octet of 
pseudo-Goldstone bosons; and / ~ /V ?a 132 MeV is the pseudo-Goldstone bosons decay 
constant. Equation (fH)| ) generates not only weak couplings of D, D* to hadronic final states 
with two pseudo-Goldstone bosons, but also the amplitudes with one light pseudoscalar 
boson in the final state; in particular, it produces the Callan-Treiman relation relating the 
form factor F ( f >7r (g 2 ) at q 2 = m? D with fr>, i.e., 

Ftijnl) ~ ^ . (18) 

At the scale we are interested in, this vertex should be corrected by a form factor that we 
call F a (q 2 ) (q is the four momentum carried by the current), 

7^-7^ 2 )- (19) 
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We do not have sufficient information on the behaviour of X a (g 2 )|, therefore, we leave it as 
a (constant) parameter and we write 

F o = 1.0±0.5. (20) 

It must be stressed, however, that, in the evaluation of the scaling behaviour (with 1/nib) of 
the charming penguin contributions, the role of this form factor is indeed relevant. Assuming 

. 3 

/ A \ 2 i 
r , as m p|, that F Q Bv (m 2 K ) scales as — ) , one gets a scaling law ml for the factorized 

contribution (in our language At+p). However, in these hypotheses, the form factor F B ' 7T (q 2 ) 

( A \ ^ 

should display a pole behaviour, to match with the — behaviour predicted by the 

\m b J 

Callan-Treiman relation at q 2 ~ m 2 B . One can now assume a similar behaviour for the form 
factor F a , which is reasonable, as the two amplitudes P — > M and P — *> MM (P heavy, M 
light mesons) derive from the same effective current ZA This implies that the contribution of 
the charming penguin contribution diagrams should be suppressed by some power 0(1/ mb) 
in the m& — > oo limit in comparison with the factorizable ones. As we are not able to define 
better the form factor F a , our evaluation should be understood as an order of magnitude 
estimate. 



3 Results 

Given these remarks, we are now ready to present our results. The Tree and Penguin 
contribution to the decay processes B — > Kw, B — > tvti, and B — > Krj^ 1 ' is obtained by the 
usual procedure of factorization using the non leptonic hamiltonian as given e.g. in ||17|| . As 



for the charming penguin contribution terms, the explicit formulae can be found in X and need 
not be reported here (in X they are denoted as Ald)- The numerical results we obtain for the 
amplitudes are reported in Table 1. We note that the phase of A^+p is due only to the weak 
interactions, while the phase in Achp is purely strong. We use the following set of parameters 
(with the notations of fl"7|l): For the Wilson coefficients [18]: c<i = 1.105, C\ = —0.228, C3 = 



0.013, c 4 = -0.029, c 5 = 0.009, c 6 = -0.033, c 7 /a = 0.005, c s /a = 0.060, c 9 /a = -1.283, 
and do/a = 0.266. Moreover, we use F S M ' (m 2 M ) w F BA/ '(0) = 0.25 (M, M' = K, vr ± )0. 
The amplitudes are evaluated using, for the CKM matrix elements, the results of the analysis 
21[: A — 0.82, p = 0.23, and 77 = 0.32. For the Kr)^ final state, we use SU(3) symmetry 



and the method of [17| with /„ = / 8 = f v = 132 MeV and 9 = 9 8 = -22°. We also notice 



tA model calculation of this form factor is in p5| . 

^This assumption is based on the dominance of the hard contribution in the QCD evaluation of the form 
factor; the actual scaling law may be affected by the behaviour in the soft region, e.g. at the end points. 
§We employ the QCD Sum Rule result of [19]; a slightly higher value is in po[. 
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Table 1: Theoretical values for A?+p (Tree+Penguin amplitude) and AchP (Charming Penguin 
amplitude). 



Process 


A r+P x 10 8 GeV 


A ChP x 10 8 GeV 


B+ -> K°7T + 


+1.69 


+2.06 + 2.36 % 


B + _> K + 7i° 


+1.21 - 0.498 z 


+1.45 + 1.67 i 


B° -> K + 7T~ 


+1.32 - 0.634 z 


+2.06 + 2.36 i 


B° -> K°7T 


-0.921 - 0.0497* 


-1.45 - 1.67i 


5+ -> 7T + 7T 


-1.35 - 1.79i 





5° -> 7T+7T- 


-1.85 - 2.16 z 


-0.576 - 0.648 z 


E° -> 7T 7T 


+0.0516 + 0.379 z 


-0.576 - 0.648 z 


-> fsT+r/ 


-0.0491 ■ 0.415 % 


+0.0830 + 0.0896 z 


5+ -> K+r/ 


+1.40 - 0.261 i 


+2.53 + 2.83 i 


5° -> fsT°r/ 


+ 0.172 - 0.0418* 


+0.0830 + 0.0896 z 


5° -> i^r/ 


+1.54 • - 0.0269i 


+2.53 + 2.83 i 



that our phase convention is such that the amplitude A{B + — > K°ir + ) differs by a sign from 
the result of p2| ; for B — > 7r 7r°, the statistical factor 1/2 in the branching ratio takes into 
account the identity of the final mesons. From the results in Table 1, we can compute the 



Table 2: Theoretical values for the CP averaged Branching Ratios (BR) compared with experi- 
mental data. Data are averages p3| from among CLEO p^j , BaBar p5| 1, Belle [^6| except for the 
upper limit that comes from [24]. 



Process 


BR XlO 6 (T+P) 


BR xlO 6 (T+P+ChP) 


BR xlO 6 (Exp.) 


B± -> K°7T ± 


~ 2.7 


18.4 ± 10.8 


17.2 ± 2.5 


B ± -> i^TT 


~ 1.6 


9.5 ± 5.5 


12.1 ± 1.7 




~ 1.9 


15.3 ± 9.9 


17.2 ± 1.5 


5° -> K°7r° 


~ 0.75 


7.4 ± 4.8 


10.3 ± 2.5 


B ± ^ ^0 


~ 4.8 


~ 4.8 


5.6 ± 1.5 


5° -> 7T+7T- 


~ 7.2 


9.7 ± 2.3 


4.4 ± 0.9 


5° -> 7T 7T 


- 0.06 


0.37 ± 0.35 


< 5.7 



Branching Ratios (BR) and the CP asymmetries for the Kn and irir final states. The CP 
averaged Branching Ratios are reported in Table 2. In the first numerical column, we report 
the results obtained by including only the Tree and Penguin contributions, i.e., At+p] in 
the second column, we give the results obtained by the full amplitude At+p + AchP] in 
the final column, we give the available data from the CLEO, Belle and BaBar experiments. 
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The errors on the branching ratios are obtained varying independently the cut-off pn in the 
range 0.5 -j- 0.7 GeV, F a in the range 0.5 -j- 1.5, and F(\p T \) = 0.065 ± 0.035, and summing 
the errors in quadrature. We have not added the errors related with the Tree and Penguin 
contribution, arising from the CKM matrix elements and from the hadronic parameters. A 
comparison between and the first and the second column shows the importance of Aqhp for 
the Kit final state, while for the ttw final states the charming penguin contribution is either 
absent (7t ± tt°) or less important (7r + 7r~). As a matter of fact, as already observed in X, using 
SU (3) symmetry one obtains for this channel 



^4chp(-B 



o 



7T + 7T 



A ChP (B° -> K+tt' 



(21) 



^cs fir 

i.e., a CKM suppression in comparison with the Kir final state. We note a general good 
agreement with the data; the only significant difference is for the 7r + 7r~ final state, which 
in our opinion should be explained by a more refined analysis of the errors in the inputs of 
the Tree and Penguin contributions (we repeat that, for the sake of simplicity, we have not 
introduced these errors in our discussion). In any event, Table 2 show that the uncertainties 
arising from the charming penguin contribution term are rather large. 

The absorptive part of AchP, which is less sensitive to theoretical uncertainties than the 
real part, provides a strong argument for a large inelastic final state interaction phase and 
for an appreciable CP violation even in the absence of the Kit and hit elastic rescattering 
phase shift. To be more quantitative, from the results in Table 1 we compute CP violating 
asymmetries for the various channels: 

- BR(B° 



A' 



A~ 



A- 



BR(B° 

br(b® 

BR(B- 



7T 1 7T 



7T 7T 



A 



00 




BR(B~ 
BR(W 



BR(B° 
BR(W 



BR(B° 

We obtain the following results: 
BR(W 



KH°) + BR(B° 



A^ 



7T 1 7T 



BR(B<~ 



7T 1 71 



0.24 ± 0.24 



(22) 



(23) 



BR(B° -> TX+TX-) + BR(B° -> tt+tt-) 

while the asymmetries for the Kir final state are reported in Fig. 2 as a function of the angle 
7 = arg(V* b ). We have not reported the asymmetry A°~ that vanishes in our approach. 
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The regions reported in these graphs correspond to a variation of the three most relevant 
parameters affecting our numerical results, i.e., the cutoff \xg G [0.5, 0.7] GeV and the form 
factors F a in eq. (p0[) and F(\p n \) G [0.03, 0.10]. We see that the variations are rather 
large, but still compatible with the CLEO [27], BaBar [P8|| , and Belle data [P9fl , which are 
as follows p3| : 



A~° = -0.096 ± 0.119, A~ + = -0.048 ± 0.068, A°~ = -0.047 ± 0.139. (24) 
The CP asymmetries we obtain are large, about 20% or more, as shown in Fig. 2. In partic- 





Figure 2: First line, CP asymmetries A (left) and A + (right); second line .A 00 . The asymmetries 
are plotted versus the angle 7 . 

ular, we find large CP asymmetries for B° — > 7r + 7r~ decays, see Eq. (p3|). The measurement 
of the weak angle a from B° — > -k + -k~ decays could still be possible once an accurate de- 
termination of the long-distance absorptive part from B — > Kit decays was obtained. Our 
results for the asymmetries are (in absolute value) compatible with Refs. |§, j30| , which also 



obtain large CP asymmetries for B — > Kir, tttt decays in phenomenological analyses of the 
charming penguin contributions. This is in contrast with the QCD-improved factorization 
model, which predicts small CP asymmetries for B — > Kn and B — > tttt decays Jl6], |23, |31 



Let us also briefly comment on the Krj and Kt]' final states. Our results for these channels 
are reported in Table 3. For the Krj' final states, one can clearly see that the charming 
penguin contribution significantly enhances the results and may be important for producing 
a large branching ratio; however, it is also clear that some relevant further contribution 
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is still missing since by no reasonable choice of the parameters can the charming penguin 
contribution alone solve the puzzle posed by experimentally very large decay fractions. We 
refer the reader to the existing literature |32|] on this subject. 

Table 3: Theoretical values of Branching Ratios (BR) for B — > Krj^'' compared with experimental 
data fro m (a) CLEO g|; (b) average between Belle @ and CLEO @. 



Process 


BR xlO 6 (T+P) 


BRxlO 6 (Ch+T+P) 


BRxlO 6 (Exp.) 


B+ _> K+T] 


~ 0.162 


0.099 ± 0.029 




B+ K+rj' 


~ 1.83 


20 ± 10 


80.0 ± 12.2 (a) 


B° K°r] 


~ 0.027 


0.058 ± 0.016 




B° K°rf 


~ 2.00 


21 ± 10 


80 ± 15 (b) 



4 Conclusions 

In conclusion, we have extended our model of the charming penguin contributions in B — > Kn 
decays to all the significant decays with two pseudoscalar mesons in the final state. Although 
the calculation presents a number of theoretical uncertainties, it clearly shows that the effect 
of the charming penguin contribution terms is overwhelming for all the B — > Kn decay modes 
while its role is less significant in the other channels. This dominance is not parametric, i.e., 
it does not contradict the dominance of the factorized amplitude in the mj — > oo limit 
discussed by several authors in the last two years [TJJ], JIB] , |3"3|| , fl3l|] , P5| . It arises from 



the CKM enhancement of the non-factorized decay chains (]T|) and their related real parts. 
The size of these charming penguin contribution terms can be estimated by an effective 
field approach, though a complete calculation is beyond the presently available theoretical 
methods. Therefore, one cannot escape the conclusion that, in spite of the proven theorems, 
the elusive non-leptonic B-decays still maintain their secrecy. 
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